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Abstract: 

We give a simple demonstration of the formula relating the glass transition 
temperature, T g , to the molar concentration x of a modifier in two types of 
glasses: binary glasses, whose composition can be denoted by X n Y m + x M p Y q , 
with X an element of Ill-rd or IV-th group (e.g. B, or Si, Ge), while 
M p Yq is an alkali oxide or chalcogenide; next, the network glasses of the type 
A x Bi_ x , e.g. Ge x Sei_ x , Si x Ttx_ x , etc. After comparison, this formula 
gives an exact expression of the parameter (3 of the modified Gibbs-Di Marzio 
equation. 
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1 Introduction 



The nature of glass transition is complex |T[ and even today remains poorly 
understood. Numerous studies |2| have been devoted to measurements and 
understanding of the glass transition temperature, T g , which is influenced 
by experimental conditions (e.g. the cooling rate of the melt). Neverthe- 
less, even for the measurements under standard conditions, there is still no 
consensus as to which structural or thermodynamical factors are responsible 
for determining T g . Kauzmann has observed that many glass- formers have 
T g ~ \T me i ting Other attempts suggested correlation with molecular or 
atomic motions with characteristic energies that might be responsible for the 
observed values of T g 0. More recently, Tanaka discovered an empirical re- 
lationship between T g and the average coordination number m per atom in 
glass : In T g ~ 1.6m + 2.3 ||. 

Nevertheless, all these descriptions become more and more complicated if 
one starts to vary the composition x of glass forming material. Typical exam- 
ples are the chalcogenide glasses (e.g. As x Se\_ x ) or binary glasses involving 
a network former (e.g. silica) and a modifier (an alkali oxide or chalcogenide, 
e.g. L12O). Some simple rules for predicting the glass transition temperature 
under standard conditions as function of the composition x and the nature of 
the atoms involved exist [[J among which Tanaka's relation J|, but there is 
still no general formulation of T g versus x which could give the precise shape 
and explain with a mathematical model why maxima in glass transition tem- 
perature can occur (as in B2O3 based systems or in Ge x Se\- x glasses at 
x=0.33 §). 

For many structural glasses, like the alkali-borate (1 — x) B2O3 + x Li 2 
glass, or the selenium-germanium network glass Ge x Sei_ x , at least at low 
concentrations of the modifier, their glass transition temperature increases 
as even small amounts of modifier are added to the melt. In many cases the 
function, T g (x), is quasi-linear only for low modifier concentrations (small 
values of x), and in some cases (e.g. the aforementioned B2O3 and Ge x Sei- x 
based glasses) displays a peak at some value of x, after which the derivative, 
dT g /dx, changes its sign and becomes negative. An elegant explanation of 
this phenomenon, based on the evaluation of average number of constraints 
per structural unit, has been proposed for the Ge x Se\_ x alloy by J.C. Phillips 

(!)■ 

We propose in this paper a simple and general model of glass formation 
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enabling us to derive the dependence of T g on x, and a very simple formula for 



its derivative at x = 0. In the refs. 101 and 111 we have derived this formula 



and applied it to the particular case of alkali-borate glass; in this paper we 
show that it is confirmed by a large set of experimental data, concerning 
mainly network glasses. 

It seems obvious that a liquid which is about to undergo the glass tran- 
sition, during which its viscosity grows exponentially, is composed of many 
clusters, various in size and shape, agglomerating as the temperature de- 
creases. Usually the process of solidification leading to the vitreous state, in- 
terpreted as a continuous random network, (the term introduced by Zachari- 
asen, fll2| ) occurs in a small temperature range, (often a function of the 
cooling rate) , between 30 and 40K around the average value called the glass 



transition temperature T„ [13|. It is known that the process is not a regular 



second order phase transition in the thermodynamical sense [IT4 



The most important elementary process of agglomeration, which is the 
creation of a new bond (usually an oxygen, sulfur or selenium bridge) between 
two clusters, can be analyzed in the ease of binary glasses. The modifier 
transforms the atoms of the glass former in a well defined manner (which 
may change when the concentration of the modifier increases); e.g., when a 
small amount of Li 2 is added to glass, the ions of lithium transform 

the threefold boron atoms of the network former into fourfold boron atoms 
[0]; when the concentration x attains about 35%, the action of the modifier 
changes, and new lithium ions begin to depolymerize the network (by creating 
non-bridging oxygens), transforming the threefold boron atoms so that they 
have only 2 free bonds left that can participate in network formation. A 
similar modification of the atoms of glass former occurs in many other glasses; 
almost always, it amounts to the change of the coordination number of the 
elementary basic structure unit of the glass former (e.g. four-, three-, two- 
, one- and zero-fold tetrahedra in IV — VI based glasses as SiS 2 — Li 2 S 
systems, corresponding to SiS2 Q [n=4..0] units ]T5|). 



2 The model 

When the modifier's concentration is not too high, we assume that the prob- 
ability of finding an "altered" atom with coordination number m' (such as 
a germanium atom in Ge x Se\- x systems) on the rim of an average cluster 
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is equal to x, and probability of finding a "regular" atom with coordination 
number m is 1 — x (such as a selenium atom in Ge x Sei- x ). Let us denote 
the non-altered ("regular") atoms of the glass former by A, and the modified 
("altered") atoms by B. 

Our aim is to evaluate the time dependence of the fluctuations of local 
concentration (i.e. local variations of the variable x), and derive the equations 
imposing the minimization of these fluctuations. This minimum should then 



correspond to a stable or meta-stable configuration |L6], [17]] ■ With only two 
types of atoms there are three elementary processes of single bond formation, 
A-A, A-B and B-B. 

The probabilities of these processes are proportional to the products of 
the relative concentrations of the atoms involved, the statistical factors which 
may be regarded as the degeneracies of the corresponding energies, and are 
proportional to the products of the corresponding valencies, and the Boltz- 
mann factors involving the energies of respective bond formation. If we look 
at the newly created pairs, their probabilities will be: 

2 

Paai = -Q-{l-xye ^ (1) 

2mm' ___B2 
Pabi = x(l -x)e *t (2) 



Q 



,/2 



Vbbi = ~q~ x e kT ( 3 ) 



where Q is the normalizing factor given by 



Q — m (1 — x) e~~*T + 2mm x(l—x) e~"^ + m x e'^r (4) 

For the moment, we exclude the possibility of a simultaneous creation of 
two bonds leading to the formation of four- and six-membered rings; such 
posibility does exist in several glasses, and can be taken into account; then the 
number of elementary processes is bigger, but the analysis remains exactly 
the same fl2], [13| . Anyhow, the network glasses we are considering here do 
not possess such an ability. Although the overall averages remain the same, 
the new pairs may create a local fluctuation in the statistics, which can be 
evaluated as: 

= -(2p BB i +Pabi) (5) 
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Denoting the average time needed to form a new bond by r, we can 
approximate the time derivative of x due to the above fluctuation as 



dx 
~dt 



1 

~ - [X 

T 



(1) 



X 



(6) 



which may be considered a good approximation as long as the time varia- 
tion of the temperature, dT/dt, can be neglected, which is true when the 
quenching rate, q = dT/dt, satisfies: 



1 dT 1 1 
1 dt 1 t 



(7) 



This inequality is particularly well satisfied for the systems with which 
we shall be dealing. For example, borate and silicate based glasses form very 
easily, and have critical cooling rates of the order of q cr u — 10~ 4 K .s _1 [ITS 



so do other network glasses we are going to investigate 19 



Now, if we want to minimize the local fluctuations, we should equate the 
above expression to 0, which amounts to finding the stationary or singular 
solutions of the differential equation (El) This leads to: 



dx If 1 







which leads to a simpler condition: 



x(l 



X) 



mil — x) (m e kT 



■me ' k T ) + m xvm e~ kT —me~ kT 



El 



(8) 



0, (9) 



There are always two singular solutions at the points x = and x = 1; there 
can exist also a third solution, given by the following expression: 



mm e k ^ 



2 

m e k T 



2 mm' e k T — m' 2 e k % — m 2 e ki 



2 -SL 

1 ^ C hT 



(10) 



However, recalling that the variable x represents the modifier concentration, 
it is physically acceptable only if < x am < 1; moreover, we want it to be 
an attractive point, which is possible only if: 



m ' gg-gi m e 2 -e 3 

— > e ^ a nd — > e kT . 



m 



m 



(11) 
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We interpret the presence of the stable solution, x am , as the manifestation 
of tendency of the system to become amorphous. On the contrary, when 
condition ( |TTD is not satisfied, the stable (attractive) solution is found at 
x = 0, which means that at the microscopic level the agglomeration tends to 
separate the two kinds of atoms, A and B. 

Note that due to the homogeneity of the expression fllPf) , only two energy 
differences are essential here: a.\ = (E 2 — Ei)/kT, and a 2 = (E 2 — E 3 )/kT. 

Various glass formers display the tendency towards the solution x am in a 
wide range of modifier concentration, also when it tends to zero [x — > 0). In 
this limit we get one condition only, which concerns only the energy difference 
E 2 — E\ : 

771 

E 2 - Ei = kT ln{—). (12) 
m 

where To is the glass transition temperature at x = 0, i.e. in the limit when 
the modifier concentration goes to zero. 

This equation summarizes the relation between the entropic and energetic 
factors that are crucial for the glass forming tendency to appear. It tells us 
that in good binary glass formers whenever m! > m (and ln(—) > 0), one 
must have E 2 > E\, and vice versa. 

This condition is what should be intuitively expected indeed when a sys- 
tem displays the tendency towards amorphisation: it behaves in a "frus- 
trated" way in the sense that the two main contributions to the probabilities 
of forming bigger clusters act in the opposite directions. Whenever the modi- 
fier raises the coordination number [m! > m), thus creating more degeneracy 
of the given energies Ei (i.e. much more possiblities of linking two entities A 
and B) and increasing the probability of agglomeration, the corresponding 

E 2 E 1 

Boltzmann factor e~~£r is smaller than for the non-modified atoms, e~kr, 
reducing the probability of agglomeration, and vice versa. 

Recalling that the stable solution corresponding to the glass-forming ten- 
dency defines an implicit function, T(x), via the relation 

E 2 E x E 3 E 2 

m (1— x) (m e kT —m e kT)-\-mx(me —me kT ) = $(a;, T) = 0, (13) 
we can easily evaluate the derivative of T with respect to x: 




with m and m' representing the average number of free valencies of the 
two species of involved atoms. 

In the limit x = the result has a particularly simple form: 



\dT- 






-l 


_dx _ 


x=0 


[ dT _ 


T=T 



(2 (m'/m) e~kr~ 

E 2 -E 1 
kT 2 



T=T 



(15) 



Inserting condition fll2|) , we obtain a general relation which can be interpreted 
as a universal law: 



dT 



dx 



Tn 



x=0 



Hi) 



(16) 



The glass transition temperature T g in binary glasses increases with the addi- 
tion of a modifier that increases the average coordination number (m' > m), 
and decreases with the addition of a modifier that decreases the average 
coordination number im! < m). 



3 Discussion 



The formula ( |16"1) can be quite easily compared with experimental data. In 
Reference ||11|| , we have demonstrated a good agreement of this formula with 
the behavior of dT g /dx at x = for the alkali-borate glass (1— x) B 2 0^ x Li 2 
and the silicate glass (1 — x) S1O2 + xCaO. Formula (13) predicted m'/m ~ 
3/2 in the first case, which leads to a positive derivative, and m'/m ~ 1/3 
giving a negative one in the second case. 

We have tried to find, whenever possible, reported glass transition temper- 
atures of glass-forming systems which were composed of a very high fraction 
of "regular" atoms (e.g. selenium, sulphur or tellurium in the case of chalco- 
genide glasses), i.e. close to 95%, in order to approach the limit value x — > 
of formula (|T6"|). We have applied this model to simple glass forming systems, 
namely chalcogenide based glasses, for which numerous experimental data are 
available. Within the ranges (close to x = 0) found in various references (see 
Table I and II), our formula agrees very well with the experimental data for 
the selenium based glasses (Table I) . By " agreement" , we mean the fact that 
the computation of the slope of T g at x = obtained with the experimen- 
tal measured values, leads to very satisfying values of m'/m. Starting from 
pure vitreous selenium whose average glass transition temperature is 316 K 
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(obtained by extrapolation from different experimental measurements |l 
we observe that all the considered glasses display an increase of the glass 
transition temperature with increasing x . Indeed, the two valenced selenium 
atoms form a network of chains with various length p0| , hence m = 2 in 
this case. The atoms of the modifier ( Ge, Si, As ) produce cross-linking 
between the 5*e-chains, thus creating new stable structural units whose m' 
is equal to 4 (Si, Ge) or 3 (As) [^l|. Obviously, equation ([TBI ) agrees only 
for the chalcogen-rich region, in particular for germanium selenide glasses. A 
similar phenomenon is observed in the Te- and S-based glasses. The average 
value of T g of vitreous tellurium, obtained by averaging and extrapolation, 
is about 343 K [|J 



whereas the T g of pure vitreous sulphur is about 245 K 



0- 



The tables below show a satisfactory agreement with our formula, and 
prove beyond any doubt that it is not a matter of coincidence. 



Compound 


(-) 

\ m Jth 


(-) 

V / exp 


Obtaii 

X 


led with 
T 9 [K] 


AT g [K] 


Reference 


Ge x Sei— x 


2.0 


2.04 


0.05 


336 


22 


s 


Si x Sei— x 


2.0 


2.04 


0.05 


336 


22 


m 


As x Se\— x 


1.5 


1.54 


0.003 


318 


2 


Hi 


Sb x Se 1 _ x 


1.5 


1.31 


0.15 


493 


177 


an 


PxSei-x 


2.5 


2.53 


0.05 


333 


17 





Table I: Different selenium based glasses. Comparison between the theo- 
retical value ofm'/m and the experimental value deduced from the slope using 
data ofTg for the lowest available concentration x. 



Compound 


(m'\ 
\ m )th 


(-) 

\ m 

\ / exp 


Obta 

X 


ined with 
T 9 [K] 


AT g [K] 


Reference 


Si x Te\_ x 


2.0 


2.11 


0.10 


389 


46 




I 


Ge x Tei_ x 


2.0 


1.97 


0.15 


419 


76 


l|27 


1 


Ga x Te X - x 


1.5 


1.45 


0.20 


528 


185 


|8 


! 


As x S\_ x 


1.5 


1.54 


0.11 


307 


63 


w 


! 


Ge x S\— x 


2.0 


1.72 


0.10 


290 


45 


w 


! 
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Table II: Different tellurium and sulphur based glasses. Comparison be- 
tween the theoretical value ofm'/m and the experimental value ofm'/m de- 
duced from the slope using data of T g for the lowest available concentration 
x. 

Nevertheless, we should mention that in certain glasses the equation (|T6|) 
does not seem to be well satined. We believe that this problem is mostly due 
to the fact that certain modifiers exhibit already a metallic character; this 
may explain why our formula can not be applied to e.g. Al x Te\_ x glasses 
[3jfl . The covalent bonding with well-defined valencies is not well suited for 



the description of the network in this case. 

A different type of problem arises in the systems which do not form glass 
below certain minimal value of x; it is obvious that our formula does not 
apply when glass can not be formed in the limit x = 0. Nevertheless, the 
formula can sometimes be extrapolated down to x = 0, when the variation 
of T g versus x represents a straight line for greater values of x. The constant 
slope allows then a comparison with flip]). In order to test the universality of 
the relation fllTf ), experimental measurements in the very low x concentration 
range remain still to be realized (as for As x Se\_ X) Table I) and also on the 
compounds which to our knowledge have never been investigated, such as 
BmSe.gs and for which we predict T g ~ 355K, and B 05 S, 95 for which we 
predict T g ~ 27hK. 



4 Relationship with the Gibbs-Di Marzio equa- 
tion 

At the beginning of this article, we have mentioned that empirical relation- 
ships between glass transition temperature and several structural or physical 
properties have been proposed during the past ten years [[| , || • Among these, 
we would like to focus our attention on the relation of equation (|16D with 
the so-called " Gibbs-Di Marzio" equation which is particularly well adapted 
for predicting T g in chalcogenide glasses. As described above (see figure 3), 
one can consider the chalcogenide glass system as a network of chains (e.g. 
selenium atoms) in which cross-linking units (such as germanium atoms) 



are inserted. Gibbs and Di Marzio [32] have developed a theory, based on 
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equilibrium principles and assuming that glass transition is a second-order 
phase transition, which relates the increase of T g to the growing presence 
of these cross-linking agents. They have applied successfully their theory to 
explain the T g data in polymers p3l p4|. An adapted theory constructed 



by Di Marzio |35j has shown that for glass systems with some chain stiff- 
ness, the glass transition temperature versus cross-linking density X could 
be expressed as: 

where T is the glass transition temperature of the initial polymeric chain 
and k a constant. 

Later on, Sreeram and co-workers || have modified this equation and have 
expressed T g in terms of the network average coordination number < r > 
which is widely used for the description of network glasses since Phillips 
has introduced the concept of < r > in his constraint theory fl36f . These 
authors have redefined for multicomponent chalcogenide glasses the cross- 
linking density X as being equal to the average coordination number of the 
cross-linked chain less the coordination number of the initial chain, i.e.: 

X = <r> - 2 (18) 

The Gibbs-Di Marzio equation may then be rewritten as: 

Tg = 1 - /?(< r > -2) (19) 

where (3 is a system depending constant, related eventually to the bond 
interchange p7| (which is responsible for the system dependent structural 



relaxation), whereas it was suggested that the constant k is universal [p5 
Sreeram et al. fitted (least-squares fit) the constant (3 to their T g measure- 
ments and obtained the value between: 

0.55 < < 0.75 (20) 

depending on the considered system and the involved atoms. 

Let us now try to relate the modified Gibbs-Di Marzio equation to for- 
mula ( |T6| ) in the pure chalcogen limit x = 0, in order to give an analytical 
expression of (3. According to Phillips [j36H , we have expressed the average 
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coordination number < r > in terms of the coordination number of the co- 
valently bonded atoms given by the 8 — N rule (where N is the number of 
the outer shell electrons), i.e. our previously defined factors m and ml: 



< r > = m(l — x) + mix 
The slope at the origin, where x = (and < r >= m) is then: 

d T g ] = T 

d < r > 



(21) 



(22) 



<r >=m (ml - m) 
which corresponds to a constant slope and to the straight line with equation: 

1 



T 



Tn 



1 + 



ml — m) In 1 



-r(< r > -m) 



(23) 



In the vicinity of the pure chalcogen region, the Gibbs-Di Marzio can be 
expressed as: 



1 9 — 



1 +0 (< r > -2) 



(24) 



which leads by identifying the two latter equations to an analytical expression 
of the constant j3, involving only the coordination number ml and m. In 
Ref. 0, Sreeram started their investigation of multicomponent chalcogenide 
systems from pure vitreous selenium, hence m = 2: 



— = (ml — m) ln( — ) 

p m 



(25) 



System 


I 3 fit 


Correlation 


Reference 




coefficient 




Ge x Se\- x 


0.74 


0.993 


i 


Ge x Sei— x 


0.72 


0.988 


I 


Ge x Sei— x 


0.65 


0.993 


PD 


Ge x Si_ x 


0.73 


0.998 


[Ml 


Si x Se,\— x 


0.81 


0.997 





Table III: Computed values of (3 obtained from a least-squares fit, for 
different binary glass systems, compared to the value (3 = 0.72 calculated 
from equation Pffij. 
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The value of (3 can now be computed for binary systems if the coordination 
number of the atoms are known, e.g. for chalcogenide based glasses, the 
possible values for (3 are 0.36 (for m' — 5), 2.47 (for m' = 3) and 0.72 (for 
m' = 4). The latter situation corresponds to the glass Ge x Se\_ x and the 
agreement of (3 — l/2/n(2) = 0.72 with the value obtained by a least-squares 
fit of the glass transition temperatures data versus < r >, is very good. 
Other IV-II systems behave very similarly, as seen in Table III. 



5 Summary and conclusion 

We have derived in this paper an analytical formula from the statistical 
consideration of the agglomeration of coordinated entities (which may be 
identified with atoms or valenced clusters), which gives the slope of the glass 
transition temperature in very poorly modified binary systems X X Y\_ X with 
x ~ 0: 



dT 



dx 



= ~^tt (26) 



The comparison of this formula with previously computed relationships of T g 
with the concentration x or the average coordination number < r >, leads to 
an exact value of the /3 parameter of the modified Gibbs-Di Marzio equation 
of glass transition temperature. 

i = (m'-m)Zn(^) (27) 

The formula seems to agree with systems investigated elsewhere. 

We believe that such an attempt should be generalized for every con- 
centration in binary glasses, in order to explain by means of the involved 
coordination numbers of the characterisitic clusters, the extrema of T g in dif- 
ferent systems such as Ge x Se\_ x or B 2 0^ based glasses. Other direction of 
investigation should be the derivation of a similar formula for pseudo-binary 
or ternary glass systems, for which there is even more experimental data 
available. This work will be the subject of forthcoming papers. 
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